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ABSTRACT: We show that for an observer in translational, rotational or gravitational 
motion, a linearly polarized plane wave has two modes of propagation in a stationary, 
homogeneous and isotropic medium according to Hertz's version of Maxwell's theory. The 
first mode is characterized by polarization at right angles to the direction of propagation 
and has a phase velocity that is controlled by the material constants of the medium. 
The second mode is characterized by polarization along the propagation direction and 
has a phase velocity that is controlled by the motion of the observer. We outline some 
applications of the second mode in emerging technologies. 

1 Introduction 

We have shown previously that in Hertz's version of Maxwell's theory, an electric field D(r, i) is 
described by the equations 

V • D = a (la) 

V x H = dD/dt + j (lb) 

j = (uV)D (lc) 
while a magnetic field B(r,i) is described by the equations^] 

V • B = 8 (2a) 



V xE = -dB/dt-b (2b) 

b = (u • V)B (2c) 

where u = dr/dt is the velocity of point r where one is observing the fields at time t and may be 
loosely described as the velocity of an observer, detector or probe at this point [[[]. Our purpose 
in this paper is to use these equations to study the propagation of light relative to accelerated 
observers (or more precisely, accelerated points of observation) in a stationary, homogeneous and 
isotropic medium with constitutive equations 

D = eE, B = (3) 

where e and /x are independent of both time t and position r. 

We are motivated to study this problem for theoretical and practical reasons. On the theo- 
retical side, we need to find experimentally testable consequences of ( |Tc| ) and (^) because these 
equations represent the point of departure of Hertz's theory from Maxwell-Einstein electrodynam- 
ics of moving bodies. Also, a rigorous classical treatment of the problem has until now seemed 
to be elusive, leaving many classical physicists to conjecture the validity of the galilean law of 
velocity composition for translating, rotating and gravitating observers (see, for example, [0). 
This conjecture can clearly not be assumed to have unlimited validity since an electromagnetic 
wave is governed by the laws of electrodynamics while a point particle for which the galilean law 
is manifestly valid is not governed by such laws. It turns out that when the limits of validity 
of this conjecture are properly studied, one finds that, unusual as it may seem, it is possible for 
light to propagate with a velocity that depends only on the motion of the observer and not on the 
permittivity or the permeability of a medium. This peculiar mode of propagation has a number 
of potential applications, some of which will be highlighted in this paper. 

Even without these novel applications, we have a practical motivation in the fact that light 
propagation relative to accelerated observers can be used to detect acceleration in devices such as 

1 Equation (|2a|) will not be used anywhere in this paper; it is included here only for reasons of completeness. 
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accelerometers (for translation), fiberoptic gyroscopes (for rotation), gravimeters and gradiometers 
(for gravitation). These devices are used in many major industries. For examples, accelerometers 
are used extensively in the automotive industry, fiberoptic gyroscopes are essential components 
of navigational systems while gravimeters and gradiometers are indispensible tools for mineral 
exploration. Hence an important goal of this work is to suggest new ways of significantly improving 
the stability and sensitivity of these devices. This goal is desirable particularly for accelerometers 
because the current technologies (capacitive, piezoelectric, piezoresistive, etc) used in these devices 
are semi-mechanical in the sense that they require the use of a proofmass whose response to 
acceleration is transduced according to Newton's second law. Although the more recent thermal 
technologies do not require a proofmass, they require a temperature-sensitive material substrate 
and therefore have many of the drawbacks of the semi-mechanical methods. By providing a useable 
closed-form solution to the problem of light propagation relative to accelerated observers, this work 
aims to provide a basis for sensing all kinds of acceleration and acceleration gradients by fiberoptic 
methods. 

We substitute (||) into ([H]) and (|t]) to get, respectively, 

VxB = ^|+j) (4a) 

VxD = -e^ + bj. (4b) 
We take the curl of (^) and use both (|8|) and ( f4a| ) to get the electric field wave equation 

c 2 [V 2 D-V(V-D)] = ^ + | + (l/ M )(Vxb), c=l/V^. (5) 

Solutions of this equation will be sought as plane waves of the form 

D = D e i(uJt+K ^ (6) 

where Dq,w and k arc independent of time t and position r. For waves of this form it is easy to 
show that 

V 2 D = (in ■ in)T> = ~k 2 T) (7a) 
V(V • D) = in{in ■ D) = —k(k ■ D) (7b) 
d 2 T)/dt 2 = (ioj ■ iuj)B = uu 2 n (7c) 



(u • V)(dD/dt) = (u • V)(iwD) = iw(u • V)D 
= (iu;)(u • i/e)D = — u>(u ■ k)D. 



(7d) 



We shall find it convenient to introduce the quantity 

tjj = u) + u • k (8a) 
which is complex if to is complex and n is real. Writing 

ip = ipi + i"ip2, lu = uji + iu)2 (8b) 
where ipi,ip2,tt)i,U2 are all real, we have from ( |8a| ) that 

V>1 = W± + U • AC, 1p 2 — U 2 (8c) 
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and from both (|8bJ) and (g) that 

D = T> ie l{uJit+R - r \ Di(t) = D e^ 2 '. (9) 

We shall characterize this plane wave by its phase velocity v p , its group velocity v g , and the time 
t s required by its amplitude to decay by a factor of 1/e (called the skin time). These quantities 
arc defined by 

v p = [u>i/k)k (10a) 

Vg=V K U) 1 (10b) 

t s = l/«a (10c) 

where k is a unit vector in the direction of k. We shall not consider the situation when u> is 
real and k is complex in this paper, but the reader is encouraged to analyze this situation as an 
exercise. 

2 Translating Observers 

We consider an observer in translational motion such that u is a function of time only, i.e., 

u = u(t), a{t)=du/dt. (11) 

Using ( |ll"| ) and (|lc|), we have 

dj/dt = d[(u • V)D]/St 

= {[du/dt] • V) D + (u • V)[dD/dt] 

= (a-V)D + (u-V)(0D/0t). (12) 

We also have 

V x b = V x {(u • V)B} by @ 

= V x {V(u • B) - u x (V x B) - (B • V) u - B x (V x u)} by (H) 
= -V x {u x (V x B)} by (g) & (0) 

= -{u(V • (V x B)) - (u • V)(V x B) + ((V x B) • V) u - (V x B)(V • u)} by @ 
= (u-V)(VxB)by @& @ 

= /i(u ■ V)(0D/0t) + M (u • V)(u ■ V)D by @. (13) 
Substituting (|l|) and (|l|) into (|) yields the wave equation 

c 2 [V 2 D - V(V ■ D)] = d 2 T>/dt 2 + 2(u ■ V)(flD/flt) + (u ■ V)(u • V)D + (a • V) D. (14) 
For the plane wave of (^|), it is easy to show that 

(a- V)D = (a-we)D 

= i(a-K)D (15a) 

(u • V)(u • V)D = (u • V)[(u • i«)D] 

= (u.«)(u-V)Dby © 
= (u ■ m)(u • i«)D 

= -(u-k) 2 D. (15b) 
Substituting ( |l5| ) and (0) into ([l4|) leads to the dispersion equation 

c [— k D + k(k ■ D)] = -u?T> - 2w(u • k)D - (u • k) 2 U + i(a • k)D (16) 
which, when simplified by taking ( |Sa| ) into account, becomes 

-c 2 k(k ■ D) = [i? - c 2 k 2 - i(a ■ k)]D. (17) 
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2.1 Transverse Waves (k. ■ D = 0) 

For these waves the left hand side of ( |l7| ) vanishes. Ignoring the trivial solution D = 0, we have 

V> 2 - cV - i(a • k) = (18) 
from which it follows that and therefore u>, is complex. We then have 

c 2 n 2 + i(a ■ k) = ip 2 



(in + # 2 ) 2 by © 



which leads to the system of equations 



with solutions 



7p 2 -ip 2 = c 2 k 2 , 2-01^2 = a- k 



= ±cnd, ip2 — ±ck\/ d 2 — 1 (a • K > 0) 



(19) 
(20) 
(21a) 



ipi = ±cnd, ip2 = d 2 — 1 (a • K < 0) 



(21b) 



where all the upper signs (or all the lower signs) are to be taken together in each set of equations, 
and d > is given by 



l + VT+tf 2 



1/2 



() 



(22) 



The corresponding frequency components as given by (^) are 

LUl = ±CKd U K, LL>2 = ±CK\/ d 2 — 1 (a • K > 0) 



(23a) 



u>i = ±cnd — U K. ll>2 = Tc K V d 2 — 1 (a- k < 0). 

Substituting these into (|l0|) gives 

v p = (±cd — u cos 9)k 
v g = ±(cdK + w)-u ^(a-/e>0) 
t s = ±(cny/d 2 - lV 1 



(23b) 



(24a) 



v p — (±cd — u cos 6)k 
v g = ±(cdd + w) — u 

t s = t(ck\/ d 2 — l) 1 

v p = (±c — u cos 6>)/5 
i> g = ±ck — u 



t s = oo 

where is the angle between u and k, and 

X. 



(a • k < 0) 



► (a-K = 0) 



w 



4cK(i 



ja-i^Kj, x = $(1 + tf 2 )~ 1/2 , = 2a«f 



(24b) 



(24c) 



a = a • k. 



(24d) 
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2.2 Nontransverse Waves (k ■ D ^ 0) 

For these waves we multiply (jlj) scalarwise by k to get 

-cV(*e • D) = [t/> 2 - c 2 k 2 - »(a ■ k)](k • D). 

Since K'D/0, we obtain 

i? = i(a- k) 

which leads to the system of equations 



with solutions 



ipf - ip 2 — 0, 2^1-02 = a- k 



tpi = ±K?y, ip2 = ±k?7 (a • k > 0) 



tpi = ±K7y, "02 = TKT] (a • K < 0) 



where 77 > satisfies 



2k 2 



The corresponding frequency components as given by (^) are 

Ll>i = ±Krj u k, u>2 — ±nr) (a • k > 0) 

Wi = ±kt; u k,, LO2 — T K V (a • ^ < 0). 
Substituting these into (|l0|) gives 



v p = (±77 — It cos 0)/? 
■u g = ±a/(4K7/) — u 
t s = ±{k V )- 1 

v p = (±r) — u cos 8)k 
v g = ±a/(4/tyj) — u 



(a - k > 0) 



> (a • k < 0) 



(25) 

(26) 

(27) 

(28a) 
(28b) 

(29) 

(30a) 
(30b) 

(31a) 



(31b) 



t s = 00 

where, as before, is the angle between u and k. 



(a • k = 0) 



3 Rotating Observers 

We consider an observer in rotational motion with velocity 

u = fJxr, n = Sl(t) 



(31c) 



(32a) 
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and accelerations 

K(t)=dSl/dt, a = au/a = Axr. 
It is easy to show from these equations that 

V • u = 0, V x u = 2fi 

and for any B, 



(32b) 



(33a) 



(B • V) u = (B • V) (□ x r) by (|32a|) 

= ((B • V) n) x r + n x ((B • V) r) 
= ft x ((B • V) r) by (|32a|) 

= ft x b. 

Taking the curl of b gives 

V x b = V x ((u ■ V)B) by © 

= V x [V(u • B) - u x (V x B) - (B • V) u - B x (V x u)] by (||) 



(33b) 



= V x ((V x B) x u) - V x ((B • V) u) - 2V x (B x ft) by (|86|) & (|33a|) 
= (V x B)(V • u) - ((V x B) • V) u + (u • V)(V x B) - u(V • (V x B)) 

- V x ((B ■ V) u) - 2V x (B x SI) by @ 
= - ((V x B) • V) u + (u ■ V)(V x B) - V x ((B ■ V) u) 

-2Vx(Bxfi) by (H|) & ©. 



(34) 



Using (33b) to rewrite the terms ((V x B) • V) u and (B • V) u respectively as fl x (V x B) and 
$1 x B, we have 

V x b = (V x B) x ft + (u ■ V)(V x B) + V x (ft x B) 
= (V x B) x ft + (u • V)(V x B) 

+ [«(V • B) + (B ■ V) ft - B(V • n) - (n • V) B] by (H) 



= (V x B) x ft + (u • V)(V x B) + n(V • B) - (ft • V) B by ( |32a|) . 



(35) 



This equation is valid for any B provided u is given by (p2|). In particular, replacing B with D 
gives 



V x j = (V x D) x ft + (u • V)(V x D) + n(V • D) - (SI ■ V) D 
while replacing B with V x B gives 

V x ((u • V)(V x B)) = (V x (V x B)) x ft + (u • V)(V x (V x B)) 

+ n(v • (v x b)) - (n • v) (v x b) 

= -Ox(Vx(Vx B)) + (u • V)(V x (V x B)) 
-(«-V) (VxB)by(H). 



(36a) 



(36b) 



Our goal is to obtain an expression from V x b in terms of V x B so that B can be eliminated 
with (Maf). For this reason we need to eliminate the last two terms of (pq). Taking the curl of both 
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sides of this equation, 

V x (V x b) = V x [(V x B) x n] + V x [(u • V)(V x B)] 

v x [n(V'B)]-vx [(n- v)B] 

= V x [(V x B) x SI] + V x [(u • V)(V x B)] 



-Six (V(V ■ B)) - (« • V) (V x B) by fl39|) & (g2aj) 

= (v x b)(v ■ si) - ((v x b) • v) si + (si • v) (v x b) - n(v • (v x 

+ V x ((u • V)(V xB))-nx (V(V • B)) - (Si • V) (V x B) by 
= V x ((u- V)(V x B)) - SI x (V(V • B)) by (§2a|) & @ 
= -SI x (V x (V x B)) + (u • V)(V x (V x B)) 



B)) 



- (SI ■ V) (V x B) - SI x (V(V • B)) by ( §6g ). 

Multiplying scalarwise by SI gives an expression 

SI ■ [V x (V x b)] = SI ■ {[(u- V)(V x (V x B))] - [(n • V) (V x B)]} 

from which B can be eliminated. To eliminate B, we use ([la|) and ( |Tc| ) to get 

(n • V) (V x B) = n (SI ■ V) [dB/dt + (u • V)D]. 

Similarly, by writing II = V x D and using (Q), 

(u • V)(V x (V x B)) = /i(u ■ V)(dU/dt) + fJ,(u • V)(V x j) 

= /i(u ■ V)[dU/dt + (u • V)n - (SI ■ V) D] 



+ /x(u • V)(n x SI) + fi(n ■ V)[«(V • D)] by 
= /u(u • V)[dIL/dt + (u • V)n - (O • V) D] 

+ /i[((u ■ v)n) xo + nx(( u - v)n)] 

+ /i[((u ■ V)«)(V ■ D) + Sl(u ■ V)(V • D)] by © & 

= /i(u ■ v) [an/at + (u • v)n - (n • v) d] 



- ^$7 x ((u • V)n) + nSl(u ■ V(V • D)) by ( |32a|) . 

It follows from @, @, @ and ® that 

n.[Vx(Vxb)]=/ifi - {Ji(u ■ Vcr) + (u • v)[A u n - (n • v)d] - (n • v) a u d} 

where we have introduced the convective or hertzian operator 

a 



A^- + (u-V) 



for convenience. 

3.1 The Wave Equation 

From (H) and dla), we have 



(37) 



(38) 



(39) 



(40) 



(41a) 



(41b) 



(l/n)(V x b) = c 2 (V 2 D - Vct) - d 2 B/dt 2 - dj/dt. 



(42) 
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Let us take the curl of this equation to get 

(l//i)(V x (V x b)) = c 2 V x (V 2 D - Va) - d 2 U/dt 2 - <9(V x j)/dt 



c 2 v x (v 2 d) - d 2 u/dt 2 - !L 



n x ft + ( u - v)n 



dt 

+ aft - (fl • V) D by (|6a]) & (|36|) 

c 2 V x (V 2 D) - d 2 U/dt 2 + flx (dTL/dt) + A x II - (a • V) II 
- (u • V) {dU/dt) - a A - fl(da/dt) + (A ■ V) D 
+ (fl-W)(dT>/dt) by (pb|). 



From (41a) and (Pt3), we finally obtain the wave equation 



c 2 n ■ [V x (V 2 D)] = ft 



d 2 n 



crA + ft^- - A x n + (a ■ V) n + 2(u • V) 



dt 



an 

~dt 



dT> 



- (A • V) D - 2 (ft ■ V) — + 0(u • Vct) + (u • V)((u • V)n) 



(u ■ V)((fi ■ V) D) - (ft ■ V) ((u • V)D) 



(43) 



(44) 



The various terms of this equation can be evaluated for the plane wave of (n) as follow. Introducing 



we have 



/3 = ft ■ (k X D), 

fl ■ [V x (V 2 D)] = ft ■ [V x ((iK) 2 T))] 
= -n 2 fl - [VxD] 
= -n 2 fl ■ [in x D] 

= -in 2 (3 



(45a) 



(45b) 



r d 2 n 



fl 



= fl 
= fl 
= fl 



dt 2 
dt 2 



(V x D) 
(in x D) 



9 2 D 



IK X 



dt 2 

fl ■ [in x (iu>) 2 T>] 

-iu 2 (3 



(45c) 



fl ■ [crA] = fl ■ [A(V • D)] 
= fl ■ [A(in ■ D)] 
= i(n-T>)(fl ■ A) 



(45d) 
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n 



dt 



-fl"»(V.D) 
= 2 -^-(iK-D) 



<9D 



at 

= ifi 2 K • (iwD) 
= -w0 2 (b-D) 

fi • [A x n] = fl • [A x (V x D)] 
= n ■ [A x (in x D)] 
= -iA ■ [ft x (k x D)] 

n-[(a-V)n]=f2-[(a-V) (VxD)] 
= «■ [(a - V) 0« x D)] 
= ifi-[KX ((a - V)D)] 
= ifl ■ [k x ((a • i«)D)] 
= fl • [(a - k)(k x D)] 
= -(a - k)0 



n 



n 



= n 



(u-V)|(VxD) 
(u'V)^xD) 



= in ■ 



(u • V) [k x 



<9D 

= • [(u • V)(k x iwD)] 
= -wfi- [k x ((u- V)D)] 
= -wfl • [k x ((u • in)T>)] 

= —iuj(u ■ k)/3 



n ■ [(A • V) D] 



n 



(45e) 



(45f) 



(45g) 



= « ■ [(A ■ i#c)D] 

:»(A-»6)(n-D) 

n • [(n • v) (iwD)] 

= i^fl • [(SI ■ V) D] 

= iun ■ [(n ■ ik)t>] 
= -uj(n ■ D)(n • k) 

fl • [fi(u • Vcr)] = f2 2 [u • V(V • D)] 
= n 2 [u- V(i« • D)] 
= i0 2 u • [V(k • D)] 

= i0 2 u • [k x (V x D) + (k ■ V) D + (D • V) k + D x (V x k)] 
= ifl 2 u ■ [k x (in x D) + (k ■ zk)D] 
= -0 2 u • [k x (k x D) + k 2 T>] 
= -0 2 (u-k)(k-D) 



(45h) 
(45i) 



(45j) 



(45k) 
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n • [(u • v)((u ■ v)n)] = n • [(u- v)((u • v)(v x d))] 

= n-[(u- V)((u- V)(kxD))] 
= iSl- [(u- V)(k x ((u- V)D))] 
= in ■ [(u ■ V)(k x ((u ■ we)D))] 
= O ■ [(u- V)((u ■ k)(k x D))] 

= SI ■ [(k x D){(u ■ V)(u ■ k)} + (u • k){(u • V)(k x D)}] 
= -j9u ■ [V(u ■ k)] - (u • k)SI • [(u ■ V)(« x D)] 

u x (V x k) + (u • V)k + (k ■ V) U + K X (V x u)] 
- (u • k)SI ■ [k X (u- V)D] 
= -/3u ■ [n x /« + 2k x fi] - (u • k)SI ■ [k X (u ■ irc)D] by (pb| ) 
= -/5[u-(kx(1) + j(u-k) 2 ] (451) 



n • [(u- v){(n • v)d}] = n • [(u ■ v){(n • m)D}] 
= z(n • «)r2 • [(u • v)D] 

= i(n • k)SI ■ [(u • zk)D] 
= -(11 • k)(u • • D) 



(45m) 



SI ■ [(SI ■ V) {(u • V)D}] = SI ■ [(SI ■ V) {(u • £k)D}] 

= ill • [(u • «){(« • V) D} + D{fi • V(u • k)}} 
= iSl ■ [(u • k){(SI ■ i«)D} + D{0 • (u x (V x k) + (u • V)k 
+ (« ■ V) u + k x (V x u))}] 



-(SI ■ k)(u ■ k)(SI ■ D) + i(Sl ■ D)f2 • [SI x k + 2k x O] by fl33bj) 
-(n-K)(u-/€)(n-D). (45n) 



Substituting (|45|) into (E3) leads to the dispersion equation 



f3[a + i(^j 2 ~c 2 k 2 )} =i^ + iQ 



where 



a = k ■ (a + SI x u) 

£ = 2(0-D)(0-k)-^ 2 (k-D) 

C = (k ■ D)(n • A) - (A • D)(0 • k). 

3.2 Coplanar Waves (/3 = 0) 

For these waves the left hand side of ( [46a| ) vanishes and on using (^), we have 
which implies that 



(46a) 



(46b) 



(47) 



(48a) 



C + =0. 



(48b) 



I f £ = by (48a), then £ = by ( 48b ) and since j3 = for coplanar waves, the dispersion equations 
(46a) and (|47j) are trivially satisfied. To obtain a nontrivial solution, we shall require that £ 7^ 0. 
With this condition, we have from (E8h that 



(49) 
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The corresponding frequency components are given by (^) as 

wi = -u- k, uj 2 = -C/£- 

Substituting these into (|l0|) gives 



which, on using (8b), leads to the system of equations 

$ " 1>l ~ - cV - (C//3) = 



(50) 



(51) 



Equation ( |51| ) shows that when £ and £ have the same signs, the wave grows while if they have 
opposite signs, the wave attenuates. For rotation with a time independent angular velocity J7, we 
have A = 0, whence £ = and the wave does not grow or attenuate. 

3.3 Noncoplanar Waves ((3^0) 

For these waves we divide both sides of ( |46a| ) by (3 and multiply by i to get 

^2_ c V-(C//3)-i(a-(£//3)V>) -0 (52) 

(53a) 
(53b) 

(54a) 
(54b) 
(55) 

(56a) 
(56b) 



with solutions 

where d > is given by 
d = l\ 



+ - a = 



i + VTTtf 2 



1/2 



C 2 K 2 f3 2 



a 



The corresponding frequency components are given by (pq) as 



LO\ = ±CK>d — U • K, Ll>2 



± 



end (3 



Substituting these into (|l0|) gives 



v p = (±cd — u cos 6) ii 
v g = ±{cdti + w) — u 

t s = 2cn(3d/{±af3 - cnd£,) 



where 



Acnd 



a + fl x u w\K > + ra7a(D x ft) + w^fl — TZ74D 



and 



ro 3 



4cKG?/3tZ7 2 



TZ72 = i C 







r 



/3(A ■ D) + 4£(f2 • D)|r, ro 4 = { 2£fl 2 + 0(SI ■ A) }T 
a\ — 2c 2 K 2 d 2 



x = tf(i + tfT 1/2 , r 



(56c) 



4d7 2 c 3 K 3 /3 2 
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Equation ( p6aj ) shows that a wave propagating in the positive direction (i.e., corresponding to the 
positive signs) with a > 0, f3 > and £ > will grow when a/3 < end!;, attenuate when a(3 > end^ 
and neither attenuate nor grow when af3 = ckg?£. Moreover, if 

e=-4 + ?fe) 1/2 (57) 



then 7 = and d = by ( |54b|) . It follows from (|56a|) that under this condition, the phase velocity 



of the wave will become independent of c and will depend only on the projection of u on k. By 



( 16t ) and (^7|) , this peculiar situation arises for an observer rotating with a time independent 
angular velocity when £ = ck(3. For waves polarized at right angles to the wave vector (k ■ D = 0), 
this condition can be written as 

u) = oj — ck, f — cos B\ cos O2 sec #3 (58) 

where 9\, 02,03 are respectively the direction angles of $7 in a righthanded cartesian coordinate 
system formed by D, k and k x D. It should be noted that if d = holds with exactitude, then 



t s = holds exactly by (56a), which means that the wave will attenuate completely before it has 
a chance to travel a sensible distance. It follows that, in practice, one should expect d = to hold 
only approximately. 

4 Gravitating Observers 

We consider an observer in gravitational motion with acceleration 

g = -qr/r 3 (59) 

where q is independent of time t and position r. The velocity u of the observer satisfies the familiar 
equations (||, pg 152) 

hx u = z + gr/r (60a) 
u x r = h (60b) 



where z and h are independent of time t and position r. By taking the vector product of (60a) 
with h and noticing that u • h = by (|60b|), we get 



u=(zxh)//i 2 + (?rxh)/r/i 2 (61) 

from which we easily obtain 

V • u = 0, Vxu = -qh/rh 2 (62a) 

V-O = 0, Vxfi = (!1xr)/r 2 = (hxg)/2)i 2 (62b) 

where, for convenience, we have introduced the angular velocity 

fl = (V x u)/2 = -qh/2rh 2 = q(r x u)/2rh 2 . (62c) 

Since u does not depend explicitly on time by (|6l|), we have 

di/dt = d[(u- V)T)]/dt 

= ([du/dt] ■ V)D + (u- V)[8D/9t] 

= (u- V)(dT>/dt). (63) 



We remark here that the presence of 7 in the denominator of the expressio n for 1? in (54t) does not involve 



a division by zero. To see this, it is sufficient to rewritejjhe expression for d in ( 54t ) by putting 7 back into the 
radicand — it was factored out of the radicand before (B7|) was contemplated. 
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It is easy to see that ( j34| ) holds for any u and any B provided fi = (V x u)/2. Let us rewrite this 
equation here for convenience, i.e., 

V x b = - ((V x B) • V) u + (u • V)(V x B) - V x ((B • V) u) - 2V x (B x fi). (64) 



Instead of (33b), we now have 



,B.r,u l BV,^ + ^]u 



rh 2 



^{{(B • V) (l/r)}(r x h) + (l/r){(B ■ V) (r x h)} 



(r x h){B-V(l/r)} + (l/r){[(B- V)r] xh + rx [(B-V)h]} 



(B-r)(rxh) Bxh 



h 2 

g_ 

h 2 



In the same way we get the following equations 



. (k ■ r)(r x fi) 
(k ■ V) u = 2<i ^ ^ -Kxn 



(65) 



(66a) 



and 



((V x B) ■ V)u = 2 j ((V X B) ; r)(r X ° } - (V x B) x fi 



(r • V) fi = (r • V) (~qh/2rh 2 ) by (|62c| ) 
= Hz/2^ 2 ){(r.V) (h/r)} 
= (- (Z /2/, 2 ){(l/r)[(r • V) h] + h[(r • V) (1/r)]} 
= (-?/2fc a ){h[r-V(l/r)]} 
= <?h/2r/i 2 



= fi by (|62c| ) 

which we shall found useful later. Substituting (|65|) into (j6^) yields 
V x b = - ((V x B) ■ V) u + (u • V)(V x B) - 2V x 



(B • r)(r x fi) 



(66b) 



(66c) 



(67) 
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The last term of this equation is expanded as 
(B • r)(r x ft 



V x 



= [V(B • r/r 2 )] x (r x n) + (B ■ r/r 2 )[V x (r x ft)] 

= -(r x ft) x [V(B • r/r 2 )] + (B ■ r/r 2 )[r(V • ft) 

- (r • V) ft + {ft ■ V) r - fl(V ■ r)] by @ 
= -(r x O) x [V(B • r/r 2 )] - fl{B • r/r 2 ) by @) & (|B6^ ) 



(rxf))x -ttV(B ■ r) + (B • r)V 



fl{B ■ r) 



"(B ■ r) + 2(B.r)[(rxfl)xr] _ (r x ft) ^ 



B r 



ft 



2r x (r x ft) \ (r x ft) 



x V(B • r) 



ft(B • r) (r x ft) 
\ 2 ' - -^r 1 x V(B ■ r) by Q62 

0(B • r) (r x ft) 



x ^ B x (V x r) 



(B • V) r + (r • V) B + r x (V x B) } by (93) 



ft{B • r) B x (r x ft) (r x ft) 



+ 



[(r-V)B]- 



(r x ft) 



r(B-O) (r x n) x [(r- V)B] (r x ft) x [r x (V x B)] 



[r x (V x B)] 

(68) 



Substituting (J68|) into (pl\) yields 



V x b = - ((V x B) • V) u + (u • V)(V x B) - 2r(B 2 °^ 

2(r x ft) x [(r • V) B] 2(r x ft) x [r x (V x B)] 



(69) 



Again, our goal is to obtain an expression from V x b that contains only V x B. We can achieve 
this by multiplying both sides scalarwise with r x ft so that the last three terms vanish and we 
get 



(r x ft) ■ (V x b) = 



r x n) • [(u • V)(V x B)] - (r x n) ■ [((V x B) • V) u] 
rxfl)- [(u- V)(V x B)] 

- 2(r x ft) ■ [r~ 2 ((V xB)-r)(rxO)-(VxB)x ft] by (|66b| ) 
r x J!) • [(u- V)(V x B)] - 2r- 2 [r • (V x B)][(r x ft) 2 ] 

h2{(r x n) • [(V x B) x ft]} 
r X ft) ■ [(u • V)(V x B)] - 2r" 2 [r • (V x B)][r 2 tt 2 - (r • ft) 2 ] 

I- 2{n 2 [r ■ (V x B)] - (r ■ ft)[ft ■ (V x B)]} by © 
r x ft) ■ [(u- V)(V x B)] + 2r~ 2 (r ■ ft) 2 [r ■ (V x B)] 

- 2(r-ft)[ft ■ (V x B)] 

r x ft) ■ [(u • V)(V x B)] by ( |S2^ ) 
/x(r x fl) • [(u- V)(OT>/0*)] 
+ M (r x n) • [(u • V)((u • V)D)] by @ & ©. 



(70) 
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4.1 The Wave Equation 

From (^) we have 

c 2 [V 2 D — V(V • D)] = d 2 D/dt 2 + dj/dt + (1/m)(V x b) 

= d 2 T)/dt 2 + (u • V)(0D/0t) + (1/m)(V x b) by 

Multiplying both sides of this equation scalarwise by r x ft yields 

c 2 (r x O) • [V 2 D - V(V • D)] = (r x n) • [d 2 T>/dt 2 ] + (r x n) • [(u • V)(8D/dt)] 

+ (i/#xn)-(vxb)] 

which, in view of (f70|), leads to the wave equation 

c 2 (r x ft) ■ [V 2 D - V(V • D)] = (r x O) ■ [<9 2 D/<9i 2 ] + 2(r x fi) ■ [(u ■ V)(dT>/dt)] 

+ (r x O)- [(u- V)((u-V)D)]. 

The last term of this equation can be evaluated for the plane wave of (||) as 

(u • V)[(u ■ V)D] = (u • V)[(u • in)T>] 

= i(u ■ «)[(u • V)D] + iD[(u ■ V)(u • k)] 
= i(u ■ «)[(u • i«)D] + zD[u • V(u ■ «)] 
= -(u • k) 2 D + iD[u • {u x (V x k) 

+ (u • V)k + (k ■ V) u + k x (V x u)}] by @ 



= -(u • k) 2 D + iD[u • {2k x O + (k ■ V) u}] by ( pcj ) 
= -(u • k) 2 D + 2iD[u • (k x n)] 
(k • r)(r x SI) 



2iT> 



— k x n 



by (56a) 



= -(u • k) 2 D + 2ir~ 2 (K • r)[u • (r x 0)]D 
= -(u • k) 2 D + 2ir~ 2 (K ■ r)[n • (u x r)]D 
= -(u • k) 2 D + 2ir" 2 (K • r)(fl • h)D by ( |S0b| ) 



= -(u • k) 2 D - iqr~ 3 (K ■ r)D by (§2cj) 
= -(u • k) 2 D + i(« • g)D by @. 

Substituting (0) and (74) into ( f73| ) leads to the dispersion equation 

[D • (r x r2)][V> 2 - cV - i(K • g)] + c 2 (k • D)[k • (r x n)] = 0. 



(71) 



(72) 



(73) 



(74) 



(75) 



If D, r and J~i are coplanar, the dispersion equation becomes [k ■ D] [k ■ (r x ft)] = which describes 
a degenerate or nonpropagating field. To obtain a nondegenerate solution, we shall require that 
D • (r x ft) 0. In this case the dispersion equation can be rewritten as 



where 



a = k ■ g, £ = (k • D 



?A 2 + c 2 (£ 2 - k 2 ) - ia = 
k • (r x ft) 



D • (r x ft) 



(k-D) 



K-(rxh) 
D • (r x h) 



(76a) 



(76b) 



4.2 Eccentric Waves (£ = re) 



For these waves, we obtain from ( |8b| ) and (76a) the system of equations 

$ -if$ = 



(77a) 
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with solutions 



tpl = ±KT], 1p 2 = ±K?7 (« • g > 0) 



1p! = ±KT), ^2 = TKV (K ■ g < 0) 



(77b) 

(78a) 
(78b) 



where, as usual, all the upper signs or all the lower signs are to be taken together in each set of 
equations, and r\ > satisfies 



2k : 



The corresponding frequency components is given by ( |8q) as 

10\ = ±KT) U K, LO2 — ±ktj (k ■ g > 0) 

W\ = ±KT/ U K, LL>2 = T K V (^ ' g < 0). 

Substituting these into (|l0|) gives 

v p = (±rj — u cos 9)k 

V g = ±g/(4KT]) U } (K g>0) 

t s = ±(ki 1 )- 1 



(78c) 



(79a) 



(79b) 



(80a) 



v p = (±?y — u cos 6)k 

Vg = ±g/(4K7/) - U 

t s = T^r 1 



(k • g < 0) 



(80b) 



J g 
t, = 00 



U COS UK 

u 



4.3 Noneccentric Waves (£ 7^ 



(K-g = 0) 



From (Sb) and (76a), we obtain for these waves the system of equations 



2^i 1P2 = a 



with solutions 



tpl = ±cdn, ■02 = ±CKy d 2 — 7 2 (k • g > 0) 



-01 = ±cdn, 02 = TcK-y/ ^ 2 — 7 2 ( K • g < 0) 



where d > satisfies 



d = j 



1 + Vi + tf 2 



1/2 



7 2 = 



C 2 



, = 



(80c) 



(81a) 
(81b) 

(82a) 
(82b) 

(83) 
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The corresponding frequency components are given by (^) as 

uj\ = ±cdn u k, UJ2 — ±cny/ d 2 — -f 2 (k • g > 0) (84a) 



uji = ±cdn u k, u>2 — TcK\/d 2 — 7 2 (k ■ g < 0). (84b) 
Substituting these into (|l0[) gives 

v p = (±cd — u cos 9)k 

v g = ±(cdK + w)-u ^ \(k- S >0) (85a) 
t. = ± (asy/tP - 7 2 ) 

u p = (±cd — u cos 0)k 
« a = ±(ct« + w)-u U«.g<0) (85b) 



< s = T (cK\/~d 



2 — ,y2 



where 



and 



X 



Acud 



g — w\K > + tn 2 (r x h) + ra7 3 D (85c) 



2cKd/3f 2 rl r(K-D) r(At-(rxh)) 
^i" x J ' n72 = ^ _ ' W3 = 2 



,n - •W^/^ (85d) 
/3 = D • (r X h), x = ^(l + ^)-^ r = ^|^. 

5 Applications 

We have shown that for an observer in translational, rotational or gravitational motion, a linearly 
polarized plane wave has two modes of propagation as summarized in Table |l|. In the regular mode, 
the phase and group velocities of the wave depend not only on wavelength and the observer's 
acceleration but also (via c as defined by (§)) on the permittivity and the permeability of a 
medium. In the irregular or heretic mode, the phase and group velocities of the wave depend 
only on wavelength and the observer's motion but not on the permittivity or permeability of 
the medium. A direct consequence of this circumstance is that phenomena (such as the Sagnac 
effect) which result in small measurable quantities due to the largeness of c can be made much 
more accurate by using heretic waves instead of regular waves. Another consequence is that heretic 
waves radiated with a given speed are observable only by an observer moving with the same speed, 
hence they are well suited for selective broadcasting. They can be used, for example, to control 
speed on a highway if a transmitter is installed on the highway and a receiver is installed in a car 
because when (and only when) the car reaches the speed limit, the receiver will pick up the heretic 
wave from the transmitter and take appropriate action. Inasmuch as for an observer moving with 
a given velocity there exists not one heretic wave with a definite frequency and wavelength but 
a complete spectrum of heretic waves with different frequencies and wavelengths, all propagating 
with the same speed as the observer, other uses of heretic waves based on their frequency and 
wavelength are conceivable. 

It is easy to see that a plane wave that is linearly polarized in its direction of propagation will 
show hereticity relative to an observer in translational, rotational or gravitational motion. For if 
D and n are collinear, we have k ■ D ^ in (|l7|) for a translating observer, (3 — in ( ff6[ ) for a 
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motion 


mode 


regular 


heretic 


translation 


transverse 
(FD 


nontransverse 
PI) 


rotation 


noncoplanar 
© 


coplanar 
dP) 


gravitation 


noneccentric 
(1) 


eccentric 
© 



Table 1: Propagation modes of a linearly polarized plane wave for an accelerated observer. The 
equation describing each mode in this paper is shown in parentheses. 

rotating observer, and £ = n in (|76| ) for a gravitating observer. In the same way, it is easy to see 
that if the wave is linearly polarized at right angles to its direction of propagation, it will behave 
in a regular way. This implies that if the angle between the polarization and the propagation 
direction is varied from ir/2 to 0, the wave speed will vary from a value controlled by the material 
constants of the medium to a value controlled by the motion of the observer. Hence if the observer 
is at rest, the wave will effectively be stopped from propagating. Propagation can be of course be 
resumed by varying the polarization angle in the opposite direction; that is, from to tt/2. 

We have shown that waves polarized at right angles to the wave vector can be heretic relative 
a rotating observer if (^) is satisfied. This result provides an interesting explanation of the 
apparently spurious signals observed in the experiment of Brillct and Hall To see this, suppose 
that D and k are on the horizontal plane so that n x D is vertical. Equation ([38]) shows that 
when the table used in the experiment is very nearly horizontal so that 9% is close to 0, both 
9\ and 02 approach tt/2. Hence cj approaches zero and no spurious signal should be observed. 
As the rotation axis of the table deviates from the vertical in such a way that Q\ ^ 7r/2 and 
#2 tt/2, #3 approaches tt/2. Hence sec #3 and co become increasingly large. This explains why 
the experimenters were able to reduce the frequency of the spurious signal from lKHz to 17Hz by 
increasingly aligning the rotation axis of the table vertically. Equation also shows that for a 
given alignment of the rotation axis, the signal should have an amplitude ~ 2f2. This means that 
the signal should peak twice for each rotation of the table, as was observed in the experiment. It 
appears therefore that Brillet and Hall may have unknowingly detected heretic waves. But it is 
not clear how the larger observed signal of 200Hz can be explained by heretic waves. For this 
and other reasons, the evidence is not as conclusive as one would like, and the spurious signals 
may very well be due to other causes [|[ p|. 

Except for the situation described by (p8|), heretic modes of propagation arise mainly because 
we have assumed a nonzero charge density in the region of interest (rr ^ 0). As a result, they may 
be less convenient than regular waves for some applications. One area where regular waves have 
found common use in industry is in the detection of acceleration. Our analysis of regular waves 
can be used to measure acceleration in three ways. First, one can use the equations for the phase 
velocity to calculate the number N of fringes obtainable from the interference of regular waves. 
Since this number depends on d which in turn depends on the component of the acceleration in 
the direction of propagation via this component of the acceleration can be calculated from 
a measurement of N. Second, one can use the equations for the group velocity of a regular 
wave to calculate the aberration ip of the wave. The aberration will be found to depend on d 
so that the component of the acceleration in the propagation direction can be calculated from a 
measurement of ip. Finally, one can measure the skin time t s directly and calculate the relevant 
acceleration component from the dependence of t s on the component. To illustrate the t s method 
for translation, we observe that to a second order approximation in ( p2| ) gives 

d 2 - 1 = i? 2 /4 
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while from (24a), we have 

d 2 - 1 = (cKt s y 2 

and from (f22|), to all orders of accuracy, 

= fl||/(c 2 K) 

where au is the component of the acceleration a parallel to k. It follows from these equations that 

2c 

a \\ = J- 

which allows au to be calculated from a measurement of t s . As is usually done in practice with 
other methods, the components of a along three orthogonal directions can be measured by using 
a system of three waves propagating in these directions. 

For an observer in uniform translational motion, it appears that one should be able to use the 



phase velocity in (24c) to calculate the number N of fringes involved in the interference of two 
regular waves, and that, by measuring N and using its dependence on the observer's speed, the 
speed can be estimated. This is indeed correct. Unfortunately, as is well known, attempts to use 
this technique with the Michelson-Morley interferometer have failed to date. This circumstance 



is interpreted in special relativity as indicating the incorrectness of the phase velocity in (24c). 
We are compelled to interprete the situation differently, however, because the phase velocity in 
question follows unavoidably from our field equations, and our field equations have been shown to 
be deductive consequences of conventions ||, |?| , which means that they cannot be falsified by any 
experiment consistent with those conventions (||, pg 214). This problem will be treated in detail 
elsewhere. We conclude with the remark that the results obtained in this paper have been shown 
to be valid for linearly polarized plane waves and may not be applied to nonplanc, unpolarized 
or nonlinear ly polarized waves without further analysis. They may also not be applied without 
further analysis in a medium that is moving, anisotropic or inhomogcneous, or in a region of space 
where the divergence theorem does not hold or charge conservation is violated. For a comparative 
study of our treatment with others, we refer the reader to the works of Phipps || , Klauber |l(]] , Goy 
and Selleri |ll| , Rizzi and Tartaglia |l2| , Lammerzahl and Haugan |l3| , Kopeikin and Schafer [[u] , 
and the references contained in these works. 



A Appendix: Vector Identities 

The following vector identities are used in this paper (see, for example, p5|). For any analytic 



vectors e, f, g, h and scalar 0, 

V x (V<£) = (86) 

V • (V x f) = (87) 
V x (V x f) = V(V • f) - V 2 f (88) 

V x (0f) = «£(V x f) - f x (V0) (89) 

(f-V)(</>g) = g(f-V0) + 0((f-V)g) (90) 

(f-V)(gxh)=gx ((f-V)h)-hx ((f-V)g) (91) 

Vx (fxg)=f(V-g)-(f-V)g+(g-V)f-g(V-f) (92) 

V(f-g)=fx (Vxg) + (f.V)g+(g-V)f+gx (Vxf) (93) 

(e x f) • (g x h) - (e • g)(f ■ h) - (e • h)(f ■ g). (94) 
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Errata 22apr2001 

• sign correction in (|^), no further consequences 

• sign correction in (|6^), minor consequences for some equations 

• constant factor correction in ( |62b| ), no further consequences 

• added (|94|) and noted its use in (|70|) 
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Errata 07may2001 



• added footnotes to clarify (J2a|) and (57) 



• minor corrections in spelling, grammar and symbols 

• changed position on the 200Hz spurious signal in the experiment of Brillet and Hall 
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